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1. Introduction 

The Hamilton-Jacobi equation 

^(t,x)+H(t,x,V<Kt,x)) = Q (1) 

plays an important role in a large variety of mathematical and physical problems. 
Apart from analytical mechanics, it appears in description of a whole range of extended 
dissipative systems featuring nonequilibrium turbulent processes, from microscales of 
condensed matter and statistical physics through mesoscale setting of free-boundary 
fluid to macroscale cosmological evolution (see, e.g., a non-exhaustive collection of 
references in [TJ). The central issue in a study of nonlinear evolution for ([rj is to 
understand, both from the mathematical and physical points of view, the behaviour of 
the system after the inevitable formation of singularities. 

A theory of weak solutions for a general Hamilton-Jacobi equation, employing 
the regularization by infinitesimal diffusion, exists since the 1970s (21 |3j 0]. In the one- 
dimensional setting this theory is essentially equivalent to the earlier theory of hyperbolic 
conservation laws in fluid mechanics, developed in the 1950s |5j El [7] . In more than one 
dimension, however, the two theories are no longer parallel. 

The theory of weak solutions for the Hamilton-Jacobi equation is closely related 
to calculus of variations, and from this point of view one can say that introduction of 
diffusion is motivated essentially by stochastic control arguments [8]. In the present 
paper we adopt a related but somewhat complementary viewpoint, in which the 
Hamilton-Jacobi equation is considered as a fluid dynamics model, and construct the 
flow of "fluid particles" inside the shock singularities of a weak solution. 

A useful example to be borne in mind when thinking about — and arguably 
the most widely known variant thereof — is the Riemann, or inviscid Burgers, equation 
du 

— + u-Vu = 0, Vxu = 0. (2) 

at 

It is obtained for the Hamiltonian H(t,x,p) = |p| 2 /2 by setting u(t,x) = V<f)(t,x). 
Note that in (J2J) it is essential that the velocity field u is curl-free, so this model is 
indeed equivalent to the Hamilton-Jacobi equation ([T]). The Riemann equation may in 
turn be considered as a limit of vanishing viscosity of the Burgers equation 
f)n ^ 

— + u»- Vti" = fiV 2 u», Vx/ = 0, (3) 

so solutions of l^j can be defined as limits of smooth solutions to as the positive 
parameter /i goes to zero. 

The Burgers equation is in fact very special: it can be exactly mapped by the Cole- 
Hopf transformation into the linear heat equation and therefore explicitly integrated [5] . 
Nonetheless the qualitative behaviour of solutions to a parabolic regularization of $1]) 
for a general convex Hamiltonian 

Jj- + H(t,x,V^)=^ 2 r (4) 
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in the limit of vanishing viscosity is similar to that for the Burgers equation. It turns 
out that as fi — > there exists a limit <f>(t, x) = lim0 /i (t, x) which is called the viscosity 
solution. Remarkably the viscosity solution can be described by a purely variational 
construction which does not use the viscous regularization at all. Below we briefly 
recall the main ideas of this variational approach. 

Assume that the Hamiltonian function H(t,x,p) is smooth and strictly convex in 
the momentum variable p, i.e., is such that for all (t,x) the graph of H(t,x,p) as a 
function of p lies above any tangent plane and contains no straight segments. This 
implies that the velocity v = V p H(t,x,p) is a one-to-one function of p. In addition, 
the Lagrangian function 

L(t, x, v) = maxp [p ■ v — H(t, x,p)] 

under the above hypotheses is smooth and strictly convex in v, although it may not be 



everywhere finite: e.g., the relativistic Hamiltonian H(t,x,p) = y 1 + \p\ 2 corresponds 



to the Lagrangian L(t,x,v) that is defined for \v\ < 1 as — yl — \v\ 2 and should be 
considered as taking value +oo elsewhere. This will not happen if in addition one 
assumes that the Hamiltonian H grows superlinearly in \p\. 

The relation between the Lagrangian and the Hamiltonian is symmetric: they are 
Legendre conjugate to one another. This relation can also be expressed in the form of 
the Young inequality: 

L(t, x, v) + H(t, x,p) > p ■ v. 

This inequality holds for all v and p and turns into equality whenever v = \7 p H(t, x,p) 
or p = W v L(t,x,v). The two functions V p H(t,x,p) and V v L(t, x, v) are inverse to 
each other; we will call them the Legendre transforms at (t, x) of p and of v. (Usually 
the term "Legendre transform" refers to the relation between the conjugate functions H 
and L; here we follow the usage adopted by A. Fathi in his works on weak KAM theory 
[5], which is more convenient in the present context.) 

Note that if H(t,x,p) = \p\ 2 /2, then L(t,x,v) = \v\ 2 /2 and the Legendre 
transform reduces to the identity v = p, blurring the distinction between velocities 
and momenta. This is another very special feature of the Burgers equation. 

Now assume that <fi(t,x) is a strong solution of the inviscid equation (P), i.e., a 
smooth function that satisfies the equation in the classical sense. For an arbitrary 
trajectory ^(t) the full time derivative of along 7 is given by 

^dF^ = ^ + V0-7 = V0-7 - H{t, 7, V0) < L(t, 7, 7), (5) 

where at the last step the Young inequality is used. This implies a bound for the 
mechanical action corresponding to the trajectory 7: 

0(*2,7(* 2 )) <0(*i,7(*i))+ 1 2 L(s,j(s),j(s))ds. (6) 



Equality in ([5]) is only achieved if 7 is the Legendre transform of V0 at (t, 7): 

7 (*) = V p tf(t,7,V0(t, 7 )). (7) 
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This represents one of Hamilton's canonical equations, with momentum given for the 
trajectory 7 by p 7 (i) := V0(t, 7(i)). The other canonical equation, p = —V X H, follows 
from (CQ) and (jZJ) because 
dVd> 

Py(t) = -Qf + (V ® V0) • 7 = -V x H{t, 7, V0) - V p # • (V <g> V0) + (V g> V0) • 7 

Therefore the bound ([6]) is achieved for trajectories satisfying Hamilton's canonical 
equations. This is a manifestation of the variational principle of the least action: 
Hamiltonian trajectories ( , y(t),p 1 (t)) are (locally) action minimizing. In particular, 
if the initial condition 

( f ) (t = 0,y) = My), (8) 
is a fixed smooth function, the identity 

0(t,z) = 0o(7(O))+ / L(s n (s),j(s))ds 



Jo 

holds for a minimizer 7 such that -y(t) = x. 

The least action principle can be used to construct the viscosity solution 
corresponding to the initial data flH]): 

<j>(t, x) = min TMt)=x (0o(7(O)) + J L{s, 7(5), 7(a)) ds) . (9) 

This is the celebrated Lax-Oleinik formula (see, e.g., [10] or (9]), which reduces a PDE 
problem (jTJ), (jSJ) to the variational problem where minimization is extended to all 
smooth curves 7 such that ~f(t) = x. 

The function defined by © is smooth at those points (t,x) where the minimizing 
trajectory is unique. In this case, the minimizer can be embedded in a smooth family 
of minimizing trajectories whose endpoints at time and t are continuously distributed 
about 7(0) and -y{t) = x. A piece of initial data 0o gets continuously deformed according 
to ([5]) along this bundle of trajectories into a piece of smooth solution to ([1]) defined 
in a neighbourhood of x at time t. Of course the Hamilton-Jacobi equation is satisfied 
in a strong sense at all points where <fi is differentiable. 

However, the crucial feature of OH]) is that generally there will be points (t, x) with 
several minimizers 7^ which start at different locations 7j(0) but bring the same value 
of action to x — 7j(t). Just as above, each of these Hamiltonian trajectories will be 
responsible for a separate smooth "piece" of solution. Thus for locations x' close to x 
the function will be represented as a pointwise minimum of smooth pieces 0,: 

0(t,£c') = minj0j(t, x'). 

As all 7j have the same terminal value of action, all these pieces intersect at (t,x): 
4>i(t,x) = (p2(t,x) = ... = <p(t,x). Thus the neighbourhood of x at time t is 
partitioned into domains where coincides with each of the smooth functions 0j and 
satisfies the Hamilton-Jacobi equation (CQ) strongly. These domains are separated 
by surfaces of various dimensions where two, or possibly three or more, pieces 0j 
intersect and hence <b is not differentiable. Such surfaces are called shocks. Note that 
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a function defined by the Lax-Oleinik formula is continuous everywhere, including 
the shocks; it is its gradient that suffers a discontinuity. In general, there are infinitely 
many continuous functions that match the initial condition (jSJ) and are differentiable 
and satisfy the Hamilton-Jacobi equation (pQ) apart from some shock surfaces, just 
as 0. A standard one-dimensional example of such nonuniqueness is provided by 
4> a (t,x) = min(a|a;| — a 2 t/2,0), which for any a > satisfies the initial condition 
4> a (0,x) = and the equation d<f) a /dt + \dcf) a /dx\ 2 /2 = apart from the shock rays 
x = ±at/2, x = 0. What distinguishes the function defined by ([9]) from all these "weak 
solutions," and grants it with important physical meaning, is that appears in the limit 
of vanishing viscosity for the regularized equation (|4j) with the initial condition (jSJ) (see, 
e.g., 0). 

For a smooth Hamiltonian it can be proved that once shocks are created they never 
disappear, although they can merge with one another. Another important physical 
feature of viscosity solutions is that minimizers can only merge with shocks but never 
leave them: all minimizers coming to some (t, x) in (Q originate at t — 0. It is easy to 
see that this is not so in the above example of <f) a , where minimizers emerge from x = 
at all times t > 0. 

Moreover, in a solution given by (EJ) a minimizer that has come to a shock 
cannot be continued any longer as a minimizing trajectory. Indeed, wherever it comes, 
there will be other trajectories originated at t — that will bring smaller values of 
action to the same location. Hence for the purpose of the least action description (j9]) 
Hamiltonian trajectories are discontinued as soon as they are absorbed by shocks. The 
set of trajectories which survive as minimizers until time t > is decreasing with t, but 
at all times it is sufficiently large to cover the whole space of final positions. 

Let us now adopt an alternative viewpoint and consider the Hamilton-Jacobi 
equation as a fluid dynamics model, assuming that Hamiltonian trajectories ([7]) are 
described by material "particles" transported by the velocity field u(t,x), which is the 
Legendre transform of the momenta field p(t, x) = V0(t, x). From this new perspective 
it is no longer natural to accept that particles annihilate once they reach a shock. 
Can therefore something be said about the dynamics of those particles that got into the 
shock, notwithstanding the fact that their trajectories cease to minimize the action? The 
problem here comes from the discontinuous nature of the velocity field u, which makes 
it impossible to construct classical solutions to the transport equation j(t) = u(t,~f). 

In dimension 1 the answer to the question above is readily available. Shocks at 
each fixed t are isolated points in the x space and as soon as a trajectory merges with 
one of them, it continues to move with the shock at all later times. This description 
is related to C. Dafermos' theory of generalized characteristics [llj which, in fact, can 
be extended to a much more general situation of nonconvex Hamiltonians and systems 
of conservation laws. However, in several space dimensions shocks become extended 
surfaces and already for equation ([1]) with a strictly convex Hamiltonian dynamics of 
trajectories inside shocks is by no means trivial. The main goal of the present paper is 
to describe a natural and canonical construction of such dynamics. 
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First results in this direction were obtained by I. Bogaevsky [T2l[T3] for the Burgers 
equation (j3j). Bogaevsky suggested to consider the transport problem for a smooth 
velocity field 

7"(*) = «"(*, V), Y(0) = y. 

Since for /i > is a smooth vector field, this equation defines a family of particle 
trajectories forming a smooth flow. The next step is to take the limit of this flow as 
\l 4- 0. Bogaevsky proved that this limit exists as a non-differentiable continuous flow, for 
which the forward derivative "y(i + 0) = lim T ^ [7(£ + t) — "f(t)]/r is defined everywhere. 
If 7(t) is located outside shocks, this derivative coincides with u(t,~/(t)). Otherwise 
there are several limit values of velocity Uj, and Bogaevsky discovered an interesting 
explicit representation for 7(^ + 0): it coincides with the center of the smallest ball that 
contains all Ui. It should be remarked that uniqueness of a limit flow in the case of 
a quadratic Hamiltonian was earlier observed by P. Cannarsa and C. Sinestrari in the 
context of propagation of singularities for the eikonal equation and differential inclusions 
PH Lemma 5.6.2]. 

The original approach in [121 H3] is based on the specific properties of the Burgers 
equation and cannot be applied in the case of general convex Hamiltonians. In 
particular, the method uses the identity of velocities and momenta, which does not hold 
in the general setting. At the same time the common wisdom says that all Hamilton- 
Jacobi equations with convex Hamiltonians must have similar properties. 

In this work we consider the above strategy, consisting in the parabolic 
regularization of the Hamilton-Jacobi equation and investigation of the vanishing 
viscosity limit for the corresponding regularized flow, in the case of general convex 
Hamiltonians. We show that such a limit exists, and derive an explicit representation 
for the forward velocity of the limit flow that extends the above result for the Burgers 
equation. Yet the mechanism powering these results in the general case is completely 
different. It is based on the fundamental uniqueness of a possible limit behaviour of 7^, 
which we discuss in detail below. 

The paper is organized as follows. In Section [2] we study the local structure of a 
viscosity solution near a singularity. We also introduce the notion of admissible velocity 
at a singularity and show that it can be determined uniquely. Moreover, the unique 
admissible velocity provides a solution to a certain convex minimization problem, which 
generalizes Bogaevsky's construction of the center of the smallest ball. In Section |3] we 
demonstrate that for any nonsmooth viscosity solution there exists a unique continuous 
nonsmooth flow of trajectories tangent to the field of admissible velocities. Section H] 
contains concluding remarks. 
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2. Local structure of viscosity solutions and admissible momenta 

Let be a viscosity solution to the Hamilton-Jacobi equation ([T]) with initial data (jSJ) . 
If there is a single minimizer coming to (t, x), then <fi is differentiable at this point and 

<f>(t + t iX + £) = <f)(t, x) + ^ r + V0 • £ + • • • 

= <f)(t, x) - H(t, x, V0) r + V0 • £ + • • • , 

where dots ••• stand for higher-orger terms. If <fi is not differentiable at (t, x), this 
means that there are several minimizers , y i such that = x, each bringing to (t, x) 
a different piece <f>i of solution. Then the Lax-Oleinik formula implies that 

4>(t + r, x + £) = mhij <j>i(t + r, x + £) 

= <f>(t, x) + min;(-ifiT + pi ■ £) H , 

where Pi := V0j(t, cc) and Hi := H(t,x,Pi). 

In the latter case neither of the expressions —HiT + Pi ■ $, provides an adequate 
linear approximation to the difference (f)(t + r, x + £) — <p(t,x), but they all majorize 
this difference up to a remainder that is linear or higher-order depending on r and 
Evidently, so does too the linear form —Ht + p ■ £ for any convex combination 

P = \Pi, H = Y] XiHi 

with Aj > 0, Aj = 1. In convex analysis these convex combinations are called 
supergradients of <fi at (i, a;) and the whole collecton of them, which is a convex 
polytope with vertices (—H^Pi), is called the super 'differential of cf> [TU [15]. We use 
R. T. Rockafellar's notation d<p(t,x) for the super differential [15J. 

To avoid a possible misunderstanding it should be noted that, although uniqueness 
of minimizer coming to (t, x) implies differentiability of (f) at t and earlier times, it does 
not imply its differentiability at any t + t > t. The following example shows how this 
may happen. The function defined for r > by 

0(t + r, x + = x) - ^r 3 - 2\i\r - + r 2 ) 3 / 2 

is a viscosity solution of the equation d<p/dt + |<90/<9x| 2 /2 = that satisfies the smooth 
initial condition (f>(t,x + £) —<j)(t,x) = — f |£| 3,/2 . For r > a shock appears at £ = 0, but 
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differentiability at r = is recovered because d<f)(t + r, x) = {— 8r 2 } x [— 4r, 4r] shrinks 
to (0, 0) as t l 0. Such points (t, x) are called preshocks pQ. 

The particular case of preshocks is an instance of a general fact: replacing gradients 
with superdifferentials allows to recover continuous differentiability, but in a weaker 
sense. Namely, suppose (t n ,x n ) converges to (t, x) and the sequence (—H n ,p n ) G 
d(j)(t n ,x n ) has a limit point (~H,p). By definition of superdifferential, 

(f)(t n + r, x n + £) - <f>(t n , x n ) < -H n t + p n ■ £ H ; (10) 

passing here to the limit and using continuity of 0, we see that (—H,p) G d<p(t,x). 
Therefore the superdifferential d(p(t x) contains all the limit points of superdifferentials 
d(p(t n ,x n ) as (t n ,x n ) converges to (t, x). 

Suppose (t, x) is a point of shock where k smooth branches <pi meet. It follows from 
the above discussion that for a particle moving from a shock point (t, x) all possible 
values of the velocity v must be such that the corresponding momentum p(v) belongs 
to the convex hull of the available momenta Pi = V<pi(t,x), 1 < % < k. However, 
one can say even more. For small positive r not all the branches <pi will contribute to 
the solution at a point (t + r, x + vt), but only those of them that correspond to a 
minimum in mhij(— Hi + Pi ■ v). Denote the corresponding set of indices 

I(v) := {1 < j < k: -Hj + Pj ■ v = min(-if i + p { ■ v)}. (11) 

i 

The set I(v) has the following physical meaning: if particle moves away from a shock 
with a given velocity v then only (pj and pj for j G I(v) are relevant. Geometrically 
one can say that the convex hull of Pj, j G I{v), is the p-projection of the face of 
the superdifferential d<j)(t, x) that looks toward an infinitesimal observer who has just 
left (t, x) with velocity v. 

This implies that any possible velocity v must satisfy the following condition. 

Admissibility Condition. A velocity v* is admissible if and only if the 
corresponding momentum belongs to the convex hull of the momenta Pj,j G I(v*). 
Namely, 

jei(v*) jei(v*) 

To make the admissibility argument rigorous it is necessary to have some control 
over the remainder term in (TTOT) . A natural function class that contains viscosity 
solutions of Hamilton-Jacobi equations and in which such control is possible is formed 
by semiconcave functions [H] . We refer a reader interested in careful proofs of this and 
other convex analytic results used in this paper to monographs [HJ US] . 

Remarkably the Admissibility Condition determines the velocity v* uniquely. 

Lemma 1 (Uniqueness) Let <fi be a viscosity solution to the Cauchy problem (QJ), 
Then at any (t, x) there exists a unique admissible velocity v* . Moreover this admissible 
velocity v* is the unique point of the global minimum for the following function 

L(v) := L(t, x, v) — mm^—Hi + Pi • v). (13) 
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Proof. Recall that L(t, x, v) is a strictly convex function of v because of assumptions 
formulated in the Introduction. Rewriting 

Li(v):=L(t,x,v)+Hi-pi-v, L(v) = max* L { (v), (14) 

we see that L(v) is a poitwise maximum of a finite number of strictly convex functions 
and therefore is strictly convex itself. Furthermore, because the Hamiltonian H(t, x,p) 
is assumed to be finite for all p, its conjugate Lagrangian L(t,x,v) grows faster than 
any linear function as \v\ increases, and all its level sets are bounded. Thus L(v) attains 
its minimum at a unique value of velocity v*. 

Let us show that this point of minimum v* satisfies the admissibility condition. 
Indeed, 

V v Li(v*) = V v L(t, x, v*) - pi = p* - pi. 

Suppose that p* does not belong to the convex hull of pj, j G I(v*). Then there exists a 
vector h such that {p*~Pj)-h < for all j G I(v*). It follows that Lj(v* + eh) < Lj(v*) 
for all j G I(v*) if e > is sufficiently small. Hence, L(v* + eh) < L(v*) for 
sufficiently small e, which contradicts our assumption that v* is a point of minimum. 
This contradiction proves that v* is admissible. 

To prove uniqueness we show that if v is admissible then it is a unique point of 
global minimum for the function L. Using the strict convexity of L, we obtain 

Lj(v + h) = L(t, x, v + h) + Hj — pj • (v + h) 

> Lj(v) + V v L(t, x,v) ■ h — pj ■ h = Lj(v) + (p — pj) ■ h, 

where p is the Legendre transform of v. Since v is admissible, p = Y^j ^jPj> where all 
Xj > and J2j = 1. Hence, J2j ^j(P~Pj)- h = KEj X j)P~J2j X jPjY h = lp-pY h = 0. 
It follows that (p — Pj) ■ h > for at least one j G I(v). Thus L(v + h) > L(v), which 
implies that v is a unique point of global minimum for L. This observation concludes 
the proof. 

The admissibility property, first formulated above in a somewhat unmanageable 
combinatorial form (fT2|) . turns out to be the optimality condition in a convex 
minimization problem given by (113p . i.e., a much simpler object. In particular, if <fr is 
differentiable at (t, x), then L(v) = L(t, x, v) + H(t, x, V</>) — V0 • v and the minimum 
in f lTHj) is achieved at the Legendre transform of V0. We thus recover Hamilton's 
equation (J7|). 

The following reformulation will clarify the connection between admissibility and 
Bogaevsky's original construction for the Burgers equation. Let Vi = V p H(t,x,Pi) 
be the velocity corresponding to the limit momentum pi and observe that Pi = 
V v L(t, x, v^. The Young inequality implies that Hi = H(t, x, p t ) = p^ ■ Vi — L(t, x, Vi) 
and therefore (fl4~l) assumes the form 

L(v) = max, [L(t, x, v) — L(t, x, Vi) — V v L(t, x, Vj) ■ (v — v^]. 

The quantity in square brackets is known as the Bregman divergence between vectors 
v and Vi, a specific measure of their separation with respect to the convex function L 
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|16j . When L(t,x,v) = \v\ 2 /2, the Bregman divergence reduces to (half) the squared 
distance between the two vectors; hence the admissible velocity v* exactly conicides 
with the centre of smallest ball containing all Vi, and Bogaevsky's result is recovered. 

Finally, we discuss the physical meaning of the function L. Consider an infinitesimal 
movement from (t, x) with velocity v. Obviously <fi(t, x) + L(t, x, v) dt — <f>(t + dt, x + 
v dt) > 0. It is easy to see that in the linear approximation in dt 

cj)(t, x) + L(t, x, v) dt - <f)(t + dt, x + v dt) = L(v) dt. 

Hence the unique admissible velocity v* minimizes the rate of the difference in action 
between the true minimizers and trajectories of particles on shocks. In other words, 
the trajectory on a shock cannot be a minimizer but it does its best to keep its surplus 
action growing as slowly as possible. 



3. The vanishing viscosity limit 

In the preceding section we constructed a canonical vector field v* = V p H(t, x,p*) 
corresponding to a given viscosity solution of the Cauchy problem (JTJ) , (JHJ) - The basis of 
this construction, the admissibility condition, appears as a natural consistency condition 
between velocities and supergradients. This condition together with the variational 
principle ( |T3l) guarantees uniqueness of the admissible pair (v*,p*). 

The vector field v*(t, x) can be decomposed into a union of smooth tangent vector 
fields defined on connected pieces of smooth shock surfaces of various dimensions as well 
as on the domain where <ft is differentiable, so dynamics in the latter domain or inside any 
piece of a smooth shock surface is given locally by a smooth flow. But globally the field 
v* is discontinuous, and it is not immediately clear if there exists an overall continuous 
flow of trajectories 7 that is compatible with v* in the sense that ^(t + 0) = v*(t,~f). 
Even less obvious is the uniqueness of such a flow. 

In order to answer these questions affirmatively we employ the vanishing viscosity 
limit for the parabolic regularization 

-^ + H(t,x,V^)=fiV 2 r, /i>0, (15) 

of the Hamilton- Jacob i equation (CQ). For sufficiently smooth initial data <fi^(t = 0, y) = 
0o (y) equation (fl5|) has a globally defined strong solution, which is locally Lipschitz with 
a constant independent of fi. Moreover, M converges as fi I to the unique viscosity 
solution <p corresponding to the same initial data. Proof of these facts may be found, 
e.g., [3], where they are established for O £ C 2,a . 
Consider now the differential equation 

Y(t) = v P H{t, y, v^(t, Y)), Y(o) = y- (16) 

For fj, > this equation has a unique solution 7^, which continuously depends on the 
initial location y. Fix a point (t, x) with t > and pick for all sufficiently small \x > 
trajectories 7 M such that 7 M (t) = x. The uniform Lipschitz property of solutions M 
implies that the curves 7^ are uniformly bounded and equicontinuous on some interval 
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containing t. Hence there exists a curve 7 and a sequence 4- sucn that lim^o 7 W — 7 
uniformly. Note that all 7 Mi and 7 are also Lipschitz with a constant independent of \x 
and that 7(t) = a;. Let furthermore v be a limit point of the "forward velocity" of the 
curve 7 at (t, x), i.e., let for some sequence r^ | 

« = lim Tfe40 — [7O + T-fc) - l{t)}- 

Of course neither the curve 7 nor the velocity v are a priori defined uniquely. 
Nevertheless it turns out that v must satisfy the admissibility condition with respect 
to the solution <p and therefore it coincides with v*. Also, trajectories of the flow 7^ 
converge as \i 1 to segments of integral trajectories of the vector field v* on smooth 
shock surfaces, establishing the uniqueness of the limit flow 7- Moreover, the limit flow 
is coalescing: if two trajectories intersect at time t, they coincide for all t' > t. All these 
statements follow from the following fact. 

Lemma 2 The flow defined by n~b}) for sufficiently small fi > collapses the 
neighbourhood of the shock surface that passes through (t,x) and is tangent to (l,v*) 
onto this surface. 

Here is a sketch of proof. Let v* be the admissible velocity corresponding to (t, x), 
p* the corresponding momentum, and define 

H* = p* ■ v* — minfjjj ■ v* — Hi). (17) 

i 

The full time derivative of the function ip^{t + r, x + £) = ^(t + r, x + £) — p* ■ £ + H*t 
along an integral trajectory 7^ of equation ffTBj) which passes through x + £ at time t + r 
is given by 

dr(t + T , r(f + T)) = ^ + 

dr at 

Convergence of viscosity solutions implies convergence of their superdifferentials (to see 
this, it is enough to replace in (fTOl) the function <p with a sequence of functions <p n 
converging pointwise). Therefore limit points of (dtp^/dt, V0 M ) belong to d<f)(t,x) for 
all (t,x), and for sufficiently small \x > 0, r > 0, £ there exists (—H fl ,p tl ) e d<p(t,x) 
and v M = V p H(t, x,p^) such that, to the linear approximation, v M = 7^ + • • • and 

^ = + H* + (p" - p*) ■ + • • • (18) 
dr 

= p'' -v* - H^ - min(p, ■ «* - if,) (19) 

i 

+ (p" - P*) • K - »*) + • • • (20) 

Lines ( TT9l) and ( )20l) contain nonnegative quantities, which are positive if p M 7^ p*. To see 
this for (1201) . observe that by the strict convexity if(£, :r,p M ) > H(t, x,p*) + (p fl —p*)-v* 
and ii(i, £C,p*) > H(t, x,p M ) + (p* — p M ) ■ and add these two inequalities. 
On the other hand, the superdifferential of the function 

lim^(t + T,X + = <f)(t + T,X + - p* ■ £ + H*T 

/40 
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at (r = 0, £ = 0) contains zero because p* corresponds to an admissible velocity. 
Therefore in the linear approximation the point £ = and other points of the shock 
surface of if) tangent to (1, v*) are local maxima of ip U P to terms of higher order, and 
the full time derivative of if) along the curve £ = v*t vanishes. As d^/dt is positive 
for trajectories that start outside the shock, we see that the flow (fToT) collapses them 
asymptotically on the shock surface. This completes the argument. 

The full details of this proof and the rigorous derivation of uniqueness of the limit 
flow 7 will be given in a forthcoming article [TTj . Here we just formulate the main result. 

Theorem 3 Let <p be a viscosity solution to the Cauchy problem for the Hamilton- 
Jacobi equation (QJ) with initial data There exists a unique flow ~f y of continuous 
trajectories such that 7 y (0) = y, ^yit + 0) is defined for all t, y and coincides with 
the admissible velocity v*(t,y y (t)) given by solution to the convex minimization problem 
for n~3\) . The trajectory -y y continuously depends on y. After a trajectory *y y comes to 
a shock, it stays inside the shock manifold for all later times. The flow is coalescing: if 
two trajectories ~f y , and "^ y „ coincide at time t, then 7 3/ /(t / ) = Tj/"(^0 f or a ^ t' > t. 

4. Concluding remarks 

Starting from a viscosity solution <f> to the Hamilton-Jacobi equation, we have 
constructed a unique continuous coalescing flow 7 y compatible with the admissible 
velocity field v* defined in Section |2] in the sense that j y (t + 0) = v*(t,j y ). This flow is 
a natural extension of the smooth flow defined by Hamilton's equation ([7]). We conclude 
with a few observations concerning this construction. 

1. Recall an important observation made in the original work of Bogaevsky: pieces 
of the shock manifold, irrespective of their dimension, are classified into restraining and 
nonrestraining depending on whether p* belongs to the interior or the boundary of 
the convex polytope formed by projection of d(f)(t, x) on the p space. Particles stay 
on restraining shocks but leave non-restraining shocks along pieces of shock manifold 
of lower codimension corresponding to faces of the boundary containing p* . Shocks 
of codimension one are always restraining; in particular, so are all shocks in the one- 
dimensional case, which makes the construction of the coalescing flow *y y trivial, as 
remarked in the introduction. Interestingly, this classification of shocks, introduced 
in [12] ("acute" and "obtuse" superdifferentials of 0), seems to have been overlooked by 
physicists before. 

2. Note that the construction of the admissible velocity v* is purely kinetic: when 
the Lagragian is "natural," i.e., has the form L(t,x,v) = K(v) — U(t,x), the value v* 
is the same for all choices of the potential term U(t, x) as long as kinetic energy K(v) 
is fixed. We owe this observation to P. Choquard. 

3. Seen as a family of continuous maps of variational origin from initial coordinates 
y to current coordinates x, the flow ~/ y is clearly relevant for optimal transportation 
problems [TH [19]. An interesting problem suggested by B. Khesin is to study the 
extremal properties of this flow. Indeed it is known from [2D] that before the first 
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shock formation the flow -y y is and action minimizing flow of diffeomorphisms, while the 
first shock formation time t* marks a conjugate point in the corresponding variational 
problem. According to the suggested view, the flow constructed above may be seen 
as a kind of saddle-point, rather than minimum, for a suitable transport optimization 
problem. 

4. Another natural context to place our construction in is that of differential 
inclusions (see, e.g., [21] )■ The flow consructed here may be seen as a solution of 
differential inclusion 

7G V p H(t,y,Pr p dcf>(t,y)), (21) 

where Pr p is the p projection of the superdifferential dip. In comparison with standard 
constructions of the theory of differential inclusions the flow -y y solves ( )2~T]) in a stronger 
sense: the forward derivative 7 y (t + 0) exists everywhere. Also, a simple modification 
of Lemma |2] gives a proof of uniqueness for inclusion ( 12TT) . 

5. The flow -y y was constructed as a limit of a parabolic regularization, and it was 
noticed above that the limit of (—d<p fl /dt,'V4> fl ) belongs to d(f)(t,x) for any (t,x). This 
statement can be refined if one considers the values of V0 M along trajectories of the 
flow ( TIT)]) . Namely, arguments of Section [3] imply that under the successive limits ji \, 
and r 4 0, the gradient (d^/dt, V0 M ) taken at (t+r, 7 m (*+t)) converges to (—H*,p*). 
Therefore 

limlim//VV(* + r,Y(t + r)) = limlim[^ + H(t, y(* + r), V0 M ) 

= H(t,x,p*)-H* 

= mm(pi ■ v* — Hi) — L(£, x, v*) 

i 

= L(v*) = min„ L(v), 

where we took into account formulas (fT7|) and fTTH]) . In other words, minimum of the 
convex minimization problem ( 1T3|) coincides with the value of the "dissipative anomaly" 
in the parabolic regularization (jl]) of the Hamilton-J acobi equation (0Q). 

6. Observe also that convergence of superdifferentials makes it possible to use 
other smoothing procedures for <p (e.g., convoluting it with a standard mollifier), 
giving the same limit "y . However, one can imagine the following completely different 
regularization of the discontinuous velocity field V p H{t,x,V 4>{t,x)). Physically 
speaking, this regularization may be characterized by a "zero Prandtl number" in 
contrast with the previous class of regularizations featuring an "infinite Prandtl 
number." 

Consider the stochastic equation 

d 7 e = V p H(t, 7 e , V0(t, Y)) dt + e dW(t), 

where W is the standard Wiener process. The corresponding stochastic flow is well 
defined in spite of the fact that V0 does not exist everywhere: whenever the trajectory 
7 e hits shocks, noise in the second term will instantaneously steer it in a random direction 
away from them. 
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Assume that as e 4- the stochastic flow 7^ tends to a limit flow 7 y , which is also 
forward differentiable. It is easy to see that, due to the averaging, the forward velocity 



where the velocities Vj(t,x) are Legendre transforms of the corresponding momenta 
Pj = V4>i(t,x) at a singular point (t,x). The coefficients 7Tj are equal to asymptotic 
values of shares of time that a trajectory 7 e spends in each of the domains where 
= (f)j. Condition ( )22|) above is another compatibility condition, in certain sense dual 
to the admissibility condition in Section |2j For want of a better term let us call such a 
velocity self- consistent. 

The self-consistent velocity is a convex combination of velocities seen by an 
infinitesimal observer leaving (t, x) with velocity . Compare this with the definition 
of admissible momentum p*, which is a convex combination of momenta see by a similar 
observer moving vith velocity v*. When H(t,x,p) = \p\ 2 /2 and v = p, self-consistent 
velocities and admissible velocities coincide. It is however clear that in the case of a 
general nonlinear Legendre transform v* ^ v* — V p H(t,x,p*). 

In view of the analogy between self-consistent velocities and admissible momenta, it 
is tempting to conjecture that the admissible velocity is also unique. Generally speaking, 
this statement is wrong, although it holds in one and two space dimensions. In higher 
dimensions there exist Hamiltonians and sets of limiting momenta for which there is 
more than one admissible velocity [17]. It is an interesting problem nevertheless to see 
whether a limiting flow still exists in the limit of weak noise in spite of nonuniqueness 
of admissible velocity. This problem carries a certain similarity with the problem of 
limit behaviour for one-dimensional Gibbs measures in the zero-temperature limit, in 
the case of nonunique ground states. 

References 

[1] Bee J and Khanin K 2007 Burgers turbulence Phys. Rep. 447 1-66 

[2] Kruzhkov S 1975 Weak solutions for Hamilton-Jacobi equations of the eikonal type Mat. Sbornik 
98 450-93 (in Russian) 

[3] Lions P-L 1982 Generalized Solutions of Hamilton-Jacobi Equations (Research Notes in 

Mathematics vol 69) (Boston: Pitman) 
[4] Crandall M G, Ishii H and Lions P-L 1992 User's guide to viscosity solutions of second order partial 

differential equations Bull. Amer, Math. Soc. 27 1-67 
[5] Hopf E 1950 The partial differential equation ut + uu x — [iu xx Commun. Pure Appl. Math. 3 201-30 
[6] Lax P 1954 Weak solutions of nonlinear hyperbolic equations and their numerical computaton 

Comm. Pure Appl. Math. 7 159-93 
[7] Oleinik O 1954 On the Cauchy problem for nonlinear equations in the class of discontinuous 

functions Dokl. Akad. Nauk SSSR 95 451-5 (in Russian) 
[8] Fleming W H and Soner H M 2005 Controlled Markov Processes and Viscosity Solutions (Stochastic 

Modelling and Applied Probability vol 25) (Berlin etc: Springer) 
[9] Fathi A 2009 Weak KAM Theorem in Lagrangian Dynamics (Cambridge: Cambridge University 

Press) 



v 




(22) 



Particle dynamics inside shocks in Hamilton- J acobi equations 



15 



[10] E W, Khanin K M, Mazel A E, Sinai Ya G 2000 Invariant measures for Burgers equation with 

stochastic forcing Ann. of Math. 151 877-960 
[11] Dafermos C M 2005 Hyperbolic Conservation Laws in Continuum Physics (Grundlehren in 

Mathematischen Wissenschaften vol 325) (Berlin etc: Springer) 



[12] Bogaevsky I 2004 Matter evolution in Burgulence Preprint math-ph/0407073 

[13] Bogaevsky I 2006 Discontinuous gradient differential equations and trajectories in the calculus of 

variations Mat. Sbornik 197 1723-51 
[14] Cannarsa P and Sinestrari C 2004 Semiconcave Functions, Hamilton- J acobi Equations, and 

Optimal Control {Progress in nonlinear differential equations and their applications vol 58) 

(Boston: Birkhauser) 

[15] Rockafellar R T 1970 Convex Analysis (Princeton: Princeton University Press) 

[16] Bregman L M 1967 The relaxation method of finding the common points of convex sets and its 

application to the solution of problems in convex programming USSR Computational Mathematics 

and Mathematical Physics 7 200-17 
[17] Khanin K and Sobolevski A 2009 On dynamics of Lagrangian particles for Hamilton- J acobi 

equations, in preparation 

[18] Gangbo W and McCann R J 1996 The geometry of optimal transportation Acta Math. 177 113-61 
[19] Villani C 2009 Optimal Transport: Old and New (Berlin etc: Springer) 

[20] Khesin B A and Misiolck G 2007 Shock waves for the Burgers equation and curvatures 
of diffcomorphism groups Analysis and singularities. Part 2, Collected papers. Dedicated to 
Academician Vladimir Igorevich Arnold on the occasion of his 10th birthday vol. 259 of Tr. 



Mat. Inst. Steklova (Moscow: Nauka) pp 77-85 (also available as preprint math/0702196) 
[21] Aubin J-P and Cellina A 1984 Differential Inclusions: Set-Valued Maps and Viability Theory 
(Berlin etc: Springer) 



